Quantum Hall Effect in Three Dimensional Layered Systems 
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Using a mapping of a layered three-dimensional system with significant inter-layer tunneling 
onto a spin-Hamiltonian, the phase diagram in the strong magnetic field limit is obtained in the 
semi-classical approximation. This phase diagram, which exhibit a metallic phase for a finite range 
of energies and magnetic fields, and the calculated associated critical exponent, v = 4/3, agree 
excellently with existing numerical calculations. The implication of this work for the quantum Hall 
effect in three dimensions is discussed. 



The quantum Hall effect is one of the hallmarks of two- 
dimensional electron systems JRe) . The possibility of the 
occurrence of the quantum Hall effect in three dimen- 
sions was explored rather early H , and precursors of the 
quantum Hall effect were observed in some three dimen- 
sional systems B . The existence of well quantized Hall 
plateaus was, however, demonstrated only in three di- 
mensional layered semiconductors with significant inter- 
layer coupling || . These layered systems have attracted 
significant theoretical interest recently, due to the pro- 
posed existence of a metallic phase for a finite range of 
energies or magnetic fields 0, and a new "crural" two- 
dimensional metallic phase on the surface W. The exis- 
tence of such a metallic phase at the surface was recently 
confirmed experimentally in measurements of the vertical 
conductance (cr 22 ) ||. 

In this work we use a mapping of the three- 
dimensional layered structure onto a two-dimensional 
spin-Hamiltonian. Using a semi-classical description we 
derive the phase-diagram Iq and obtain the critical ex- 
ponent v, describing the divergence of the localization 
length £, as one approaches the transition from the in- 
sulating side, £ ~ \E — E c \~", or £ ~ \B - B c \~ u , where 
E c and B c are the critical energy and magnetic field, re- 
spectively. The derived critical exponent v = 4/3 agrees 
excellently with existing numerical data, v = 1.35 ±0.15, 
obtained both for a layered system and three-dimensional 
tight-binding model H, and v = 1.45 ±0.25, obtained in 
Ref. || from a layered network model |1(J . 
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Fig. 1. The system studied in this work. 

We start with the Hamiltonian describing a spinless 
electron in a system of N coupled two-dimensional layers 



(see Fig. 1) 
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[(ft - eA,/c) 2 /2m + V t (x, y) + T^x, y)] 



(1) 



where x and y are coordinates in the plane and the sum- 
mation is over the layers. The first and second terms in 
the brackets describe the kinetic and potential energies 
within a layer, while the third term describe the hopping 
between adjacent layers, which may depend on the posi- 
tion in the plane. The layer potentials are assumed to be 
independently distributed with zero mean. 

We now associate with the electron a spin-index that 
corresponds to the layer index in (PP) . The interlayer tun- 
neling will now correspond to spin-raising and lowering 
operations. In order to describe the different potential 
landscape in each layer, we add a random S z term to the 
Hamiltonian that now describes a spin-,!? (= (N — l)/2) 
electron moving in two dimensions, 

H = (p - eA/cf/2m + U{x, y) + AU(x, y)S z (2) 
+ t{x,y)S+ +t*(x,y)S-. 

The second term describes a random potential indepen- 
dent of the spin (layer) index. The third term accounts 
for the different potentials for the different spin direction, 
by a random shift of the potential between adjacent lay- 
ers (at each point of the plane). Thus at each point the 
electron sees a different potential in each layer (or for 
each spin direction). Since the shift AU(x,y) is random 
in sign and in magnitude, the average potential in each 
layer is the same pd| . 

The Hamiltonian (Eh can now be simply written as 



1 



H = (p - eA/c)72m + U(x, y) + -^S • H(x, y), 



(3) 



namely a spin-S" electron moving in two-dimensions un- 
der the influence of a random potential and a random 
magnetic field (coupled to its spin). The advantage 
of this representation is that one can try to generalize 
methods that worked for the two-dimensional case, in 
the absence of a random field, to include the effects of 
the field. In the following we will concentrate on the 
large (uniform) magnetic field limit, where the kinetic 
energy is quenched and one may treat the electrons semi- 
classically. In the absence of the random field the elec- 
tron moves along equi-potential lines. As is well known 
in this case |L2| , electrons with too small an energy will 
be trapped around potential valleys, while for too high an 



energy they will be trapped around potential hills. There 
is a single "critical" energy where the electron trajectory 
percolates through the system. This corresponds to the 
quantum Hall transition, where there is a single energy 
(at the center of the Landau level in case of symmet- 
rically distributed random potentials) where states are 
extended. 

In the present case, in the same strong magnetic field 
limit, it is the total energy - the potential energy plus 
the spin energy (due to the random field) that is con- 
served. Thus, as the electron rotates its spin along the 
trajectory, it exchanges energy between the potential en- 
ergy and the spin-energy, such that the total is conserved. 
The range of potential energies accessible by the electron 
has a width A = 2Hr, where Hr is the typical amplitude 
of the random field. Consequently, even if the electron 
does not have the correct (critical) potential energy to 
percolate through the system to begin with, it can still 
do that as long as its total energy is within Hr of the 
critical energy. 




Fig. 2. The classical trajectories of an electron in a strong 
magnetic field. Without random field, the electron follows 
equipotential lines (solid curve); with increasing random field 
it explores larger portion of the potential energy landscape 
(circles), until for large enough random field, it can go through 
the saddle point (broken line). 

An example is depicted in Fig. 2. The classical equa- 
tions of motion for the Hamiltonian (0) with U(x, y) cor- 
responding to two impurities (the equipotential lines ap- 
pear as thin solid curves) were integrated. In the absence 
of a random field (a solid thick curve), the electron follows 
a single equipotential line, with superimposed cyclotron 
oscillations, and is trapped around one impurity. With 
increasing random field the electron explores a larger por- 
tion of the potential energy landscape (see, e.g., the tra- 
jectory denoted by circles), until, for large enough ran- 
dom field (broken line), the electron can go through the 
saddle point and percolate away. In the original layered 
system, this process corresponds to the possibility of the 
electron tunneling to a different layer and drifting along a 



different potential line (with the same potential energy). 
Thus as the energy is increased, before percolation occurs 
in a single layer, there will be a percolating path consist- 
ing of equipotential lines in different layers, connected by 
inter-layer tunneling events. 

Since the random magnetic field amplitude Hr ~ 
^/t 2 + (A[/) 2 , one expects a region of extended states 
that increases with t, leading to the phase diagram de- 
picted in Fig. 3. For any finite t there exists a finite 
range of energies (or magnetic fields) where the system 
is metallic. Accordingly, even at T = the transition be- 
tween Hall plateaus will not be sharp, but rather occur 
in a finite range of magnetic fields or gate voltages. 

Interestingly, in the present semi-classical description 
such a metallic phase will occur even for an infinitesi- 
mal tunneling matrix element t. The reason is that once 
t =/= the electron can, in principle, rotates its spin (tun- 
nel between layers) and explore the whole energy range 
allowed by conservation of total energy. We know, how- 
ever, that quantum mechanically, for small enough tun- 
neling matrix element, the electron will be localized in 
spin-space and the range of potential energies available 
(i.e. the width of the metallic region in phase space) will 
be much smaller than one expects classically, going to 
zero as t — » Jl3[ |. Thus, there is a region in the phase 
diagram (the shaded part of Fig. 3), where the electron 
is localized quantum mechanically, but its classical tra- 
jectory is extended. The derived phase diagram (Fig. 3) 
agrees with the phase diagram established numerically 
by Chalker and Dohmen M. 




energy 

Fig. 3. Phase diagram of the layered system. For finite inter- 
layer tunneling t there is a finite range of energies where ex- 
tended states exist. The shaded region is where the wavefunc- 
tions are localized, but the classical trajectories are extended 
(see text). 

We now turn to the critical behavior. For the 
two-dimensional quantum Hall problem Mil'nikov and 
Sokolov ITjJIq] used the following argument to predict 
the critical exponent. In the classical description, away 
from the critical energy E c , the electron is confined to 
a percolation cluster of typical size £ p , the percolation 
coherence length. Near the threshold £ p ~ \E C — E\~ Vp , 
where v p — 4/3 is the two-dimensional percolation ex- 
ponent. As one approaches the transition the clusters 



approach each other near saddle points of the potential 
energy landscape. While classically the electron cannot 
move from one cluster to another, quantum mechani- 
cally it can tunnel through the potential barrier. If the 
electron energy E is close enough to the transition, the 
potential barrier is close to parabolic and the tunneling 
probability through such as saddle point is proportional 
to Exp[— (E c — E)]. The number of such saddle points 
through which tunneling occurs in a system of length L 
is typically L/£ p . Since the transmission coefficient is 
multiplicative, the conductance (or the tunneling proba- 
bility) through the whole system is 



C2D 
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L/Z P 



= g-i/«2D 



(4) 



with &d ~ (E c - E)- U2D and v 2 d = v P + 1 = 7/3. 
The best numerical estimate of the critical exponent 
V2D = 2.35 ± 0.02 |lq ], which is supported by experi- 
mental data |L7|| , has a surprisingly excellent agreement 
with the result of the above argument, especially in view 
of the crudeness of the argument. 

This argument can be generalized to the present prob- 
lem [la ], as it is also expressed in terms of a two- 
dimensional Hamiltonian. In the presence of inter-layer 
tunneling (random field), the only difference between the 
present problem and the two-dimensional problem is the 
fact that the critical energy E c is not equal to the poten- 
tial energy of the saddle-point, but is Hr away from it. 
Thus 



vs D ~[e- H »] L/ *' = e- L 'e>°, 



(5) 



mation p9| . In this approximation one neglects the ad- 
ditional IPVU term in the parentheses, and the Hamil- 
tonian can be trivially diagonalized in spin-space. The 
random field serves as an additional potential energy, 
which is different for each spin-direction (and its average 
is proportional to HrS z ). Consequently, in this approx- 
imation one expects N separate transitions, each of the 
two-dimensional type (see Fig. 4). (Note that these tran- 
sitions are not related to the different layers, but rather to 
different coherent superpositions of the wave- functions in 
different layers). Since the separate transitions can only 
be resolved for energies smaller than A/N, one expect 
in this case a crossover from a three-dimensional critical 
behavior, for \E — E c \ > A/N to a two-dimensional crit- 
ical behavior for \E — E c \ < A/N (the two-dimensional 
behavior can only be seen for temperatures smaller than 
A/N), 



with ^3£> ~ (_E C — E)~ v and v = v v = 4/3. One finds 
the surprising result that the critical exponent for the 
quantum three-dimensional problem is equal to the two- 
dimensional classical percolation exponent. This result 
is in excellent agreement with existing numerical esti- 
mates, v — 1.35 ± 0.15, obtained both for a layered sys- 
tem and three-dimensional tight-binding model [Bl, and 
v = 1.45 ± 0.25, obtained |(| from a layered network 
model g(J. 

Consider now the Hall conductance a xy . If the inter- 
layer tunneling t is equal to zero, the system is a collec- 
tion of N independent two-dimensional layers, all with 
the same critical energy. Thus a xy will jump by e 2 jh in 
all layers simultaneously (see Fig. 4), i.e., it will have a 
single step of height Ne 2 /h (which corresponds to a con- 
ductance per layer or conductivity of e 2 /h). For finite 
t (or finite random field) the situation is quite different. 
To see this we first carry out a local SU(N) gauge trans- 
formation in spin space, to rotate the spin by a unitary 
matrix U(x,y), such that the z-direction always lies in 
the direction of the random field. This exact transfor- 
mation maps the Hamiltonian (ph onto the equivalent 
Hamiltonian Jl8| 

H = (p - eA/c - iW^U) 2 /2m + U(x, y) + ^S z \H(x, y)\ 

S (6) 

If the potential energy and the inter-layer tunneling vary 
slowly in space, one may apply the adiabatic approxi- 



i = A x e- 



7(e) 



m 



1 e<l 

A 2 e^-^n e > i 



(■7) 



with e = (E — E C )/(A/N). Thus, the effective exponent 
v will crossover from its three-dimensional (~ 4/3) to 
the two-dimensional value (~ 7/3), as one gets closer to 
the critical point from the insulating side. Interestingly, 
for the case v-id = 7/3 and v-^d = 4/3 the scaling func- 
tion /(e) may be analytic. This crossover can be studied 
via the critical behavior of the conductance (Eqs.(Q) and 
(|)), or by that of da xy /dB |(J. 
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Fig. 4. The change in the Hall conductance at the transition. 
For t = there is a two-dimensional behavior (a single step). 
For ( / 0, then in the adiabatic approximation one find a 
series of smaller steps, of the number of layers (see text). In 
the three dimensional limit, as the number of layers increases, 
one expects a smooth transition between the quantized values 
(a metallic region). 

In the adiabatic approximation there is a zero- 
temperature metallic phase only in the true three- 
dimensional limit (N — > oo), which is the classical limit 



(S — > oo) of the spin-problem (Fig. 4). As non- 
adiabaticity (the additional term in the parentheses in 
Eq.(g)) is switched on, the different spin-states that were 
the eigenstates of the system in the adiabatic limit get 
coupled. It is not clear if this coupling will smear out the 
separate transitions even for a finite number of layers. It 
is known that there may occur transitions between the 
expected adiabatic behavior to a different behavior (as a 
function of e.g. the tunneling matrix element), even for 
the two-layer problem |21| , and it remains to be seen if 
such a deviation from the adiabatic limit will also occur 
for a finite number of layers. We hope that this work will 
motivate further studies in this direction. 

To conclude, we have used a mapping onto a two- 
dimensional spin-Hamiltonian to describe the physics of 
the quantum Hall effect in three-dimensional layered sys- 
tems. This mapping was used mainly for conceptual 
reasons, in order to allow us to extend methods ap- 
plied in the traditional two-dimensional quantum Hall 
systems to the present case. The arguments presented 
here, however, could be directly applied to the origi- 
nal three-dimensional system, and thus none of the re- 
sults of this paper depends on the particular form of the 
spin-Hamiltonian. For example, in the three-dimensional 
layered system, the potential and the hopping part of 
the Hamiltonian (Eq.(Q)) can be recast in a form of a 
position-dependent N x N matrix. Diagonalizing this 
matrix locally and carrying out a unitary local rotation 
in layer-space, will lead to a Hamiltonian of the form m), 
and to all the results of the last section. Similar argu- 
ments can be made to derive the phase-diagram and the 
critical exponent. 

The author thanks A. Stern for several discussions. 
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